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Abstract 


{Q> 


= source-intensities vector 


A possible reason is suggested for the induced 
rolling moments occurring on wraparound-fin config- 
urations in subsonic flight at zero angle of attack. 
The subsonic potential flow over the configuration 
at zero incidence is solved numerically. The body 
is simulated by a distribution of sources along its 
axis, and the fins are described by a vortex- 
lattice method. It is shown that Vdiii.g moments 
can be induced on the antisymmetric i Ju . by the 
radial flow generated at the base of the configura- 
tion, either over the converging separated wake, or 
over the diverging plume of a rocket motor. 
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Nomenclature 

{x c > 

= rolling-moment coefficient, normalized by 

body diameter and body cross-section area {Xp} 


radius vector, Eq. (8) 

local body radius 

velocities in system (x,r,0) 

vector of normal velocities due to the 
uniform field U^ and body-induced flow 
field 

undisturbed uniform field velocity 

vector of normal velocities due to the 
vortices system 

sources position vector 

collocation points position vector 




P 


= roll-damp ing-moment coefficient, sec/rad 


(x,r,0) = cylindrical coordinate system 


=* rolling-moment coefficient slope due to x S p 

fin deflection, 1/deg 


= wake length up to the rear stagnation 
point 


dt 

dV 

D 


« pressure coefficient 

(x,y,z) 

= Cartesian coordinate system 

= line-vortex segment vector 

3 

*» compressibility factor 3 2 = 1 - 

= velocity induced by line vortex segment, 
Eq. (8) 

(y) 

- vector of vortex strengths, per unit 
length 

- body diameter 

f 

* vortex strength vector 


[H] 

I 


J 


= influence coefficient matrix, Eq, (9) 6 

= number of source or vortex and column 0 D 

number in the matrix of Eqs. (7) and (9) 

A LE 

= number of collocation point and line 

number in Eqs. (7) and (9) A-pj? 

= total number of sources used to describe <J> 
the body 

<J>R 

- free-stream Mach number 

= unit vector normal to aerodynamic surfaces ^ 


= fin deflection angle* deg 

=* fin opening angle, deg 

“ fin leading-edge sweep angle, deg 

- fin trailing-edge sweep angle, deg 

- potential function 

- roll angle position of sting-mounted 
model, deg 

- fin span angle, deg 


N v = total number of fin panels 

Pj/P^ = rocket motor pressure ratio 


'‘National Research Council — Research Associ- 
ate, Aerodynamics Research Branch; Associate Pro- 
fessor, Technion-IIT . Member AIAA. 

^'Graduate Student, Department of Aeronautical 
Engineering . 


Introduction 

The last decade has seen the development of 
many tube-launched missiles, from ground-launched, 
saturation field-artillery rockets to air- 
launched weapons, and from simple unguided rockets 
to sophisticated cruise missiles. Tube launching 
was chosen for its packaging convenience and to 
increase the reliability of the rocket motor. 
However, such a launching procedure requires fold- 
ing aerodynamic stabilizers because of its packag- 
ing constraints. The stabilizers would be folded 
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in the stowed position to fit within a circular 
cylinder and would deploy instantly after launch. 
Wraparound-fin (WAF) configurations (see Fig. 1) 
meet this requirement and at the same time maximize 
the volume available for the missiles subsystems, 
especially for the nozzle exit* 

It has been repeatedly observed that WAF con- 
figurations have conventional longitudinal aerody- 
namic characteristics 1- ** equal to those of config- 
urations with planar fins of identical planform. 
However, as early as I960, Feather stone 1 pointed 
out that WAF configurations displayed a seemingly 
unpredictable lateral behavior. Of special inter- 
est for the present investigation were the induced 
rolling moments. While these moments were under- 
standable (if not predictable) at the higher angles 
of attack as a result of the asymmetry of the 
curved fins, 4,5 the moments measured at zero angle 
of attack s presented an enigma. An even greater 
enigma was the roll reversal that occurred when WAF 
configurations accelerated or decelerated through 
the sonic-speed range. 6 

Because of the significance of WAF configura- 
tions in weapons technology and their unusual lat- 
eral characteristics, an international research 
program [the Technical Cooperation Program (TTCP) 7 ] 
was set up in 1969 at a meeting at Eglin AFB. In 
the ensuing research effort (encompassing the 
U.S. Armed Forces, the United Kingdom, Canada, and 
Australia) an intensive experimental investigation 
program was launched to study the aerodynamic char- 
acteristics of a standard configuration (Fig. 1) 
with 20 different fins and two different bodies, 
at Mach numbers from 0.3 to 3.0 (Refs. 8-10). 

The experimental results of this program proved 
that the static longitudinal aerodynamic character- 
istics of WAF configurations did not differ from 
those of conventional configurations with identical 
bodies and flat fins, if the flat fin area was 
equal to the projected area of the curved fins. 

This conclusion meant that WAF configurations could 
be used where flat fins had been used before, with- 
out any penalty in their longitudinal aerodynamics 
and even without a severe increase in their total 
drag. However, those results also substantiated 
Featberstone 1 s 1 findings of unusual lateral aero- 
dynamics. In spite of the relatively large scatter 
in the experimental data obtained from the sting- 
mounted models in the various wind tunnels 
involved, 10 there was no doubt that WAF configura- 
tions developed induced rolling moments, even at 
zero angle of attack. Generally speaking, these 
moments acted upon the fins in the direction of 
their center of curvature at subsonic speeds and 
changed direction in transition to supersonic 
speeds, or vice versa. This behavior could be 
detrimental to the stability of WAF missiles 
because roll resonance, or pitch-roll lock-in and 
catastrophic yaw, due to roll-yaw coupling could 
develop at the point of roll reversal, increasing 
the scatter of such missiles around their targets. 7 
Although the problem can be overcome by intention- 
ally spinning the missile (by means of canted fins, 
for example) the TTCP participants still felt the 
need to identify the origin and to understand the 
nature of the induced rolling moments. However, 
in spite of years of investigative effort, the 
answer still eluded them. 

Intrigued by this problem, the authors tried 
an analytic approach * A lifting-line theory for 


antisymmetric wings with spanwise curvature in 
uniform flow was developed, with the intent that 
the asymmetric tip vortices generated at any finite 
angle of attack would retain some vorticity when 
the angle was gradually decreased to zero. How- 
ever, nothing but the trivial solution could be 
found at zero incidence. Next, the authors modi- 
fied the vortex-lattice code for wings alone 11 to 
handle antisymmetric, spanwise-curved wings. 

Again, wraparound fins at zero incidence in uniform 
flow failed to produce rolling moments. With the 
failure of both these approaches in uniform oncom- 
ing flow, and with some experimental indication 6 
that rocket-motor burnout also gave rise to roll 
reversal, the authors suspected that, in subsonic 
flow, the rolling moments could be generated by a 
radial-velocity component induced on the fins by 
the wake flow field. 

The purpose of this work was to investigate 
this idea and to prove its feasibility. It was 
not meant to be a rigorous solution of the flow 
field around a WAF configuration at any given angle 
of attack. Consequently, the authors decided to 
use existing approximate numerical methods with 
some simplifying assumptions in order to obtain at 
least a qualitative answer at a minimum effort and 
cost. 


The Mathematical Model 

The mathematical model, as well as the numeri- 
cal methods used here, is well known and is 
described only briefly. The interested reader 
may consult the relevant references for more 
details . 

Since there is no reason to expect viscous 
effects in the mechanism generating the rolling 
moments, it is sufficient to describe the potential 
flow about the WAF configuration. If the flow is 
incompressible it can be described by the well- 
known Laplace equation for the velocity potent! i 

V 2 <j> = 0 (1) 

where is the potential function, and by the 
tangency boundary condition on all solid surfaces 
of the configuration: 


where n is the unit vector normal to the surface* 
Further boundary conditions can be imposed such as 
the Kutta condition at the fins' trailing edge, 
and/or a prescribed wake geometry (linear wake), or 
a force-free, rolled-up wake (nonlinear wake). 

Such a problem can be solved in a given uniform 
flow by an assumed distribution of a discrete 
number of potential elemental functions (i.e., 
sources, doublets, and vortices) inside the body 
or on its surfaces. The intensities of these 
functions are determined by satisfying the bound- 
ary conditions. With the intensities known, the 
induced perturbation velocities at every point in 
the flow field and the pressure distribution over 
the configuration's surfaces can then be computed. 

If the solution has to be extended into the 
Subsonic compressible flow regime, as is done in 
this work, then the linearized, small-perturbation 
potential equation 
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that approximately describes the compressible flow 
is transformed into the incompressible potential 
Eq. (1) by the Gothert transformation: 

x « x ? /i3 , y *= y* , z * z' (4) 

where the primed quantities belong to the physical 
compressible flow and the unprimed ones to the 
transformed incompressible flow, and the same solu- 
tion procedure, described above, is used. With 
this approach, what remains to be decided is the 
exact model of the potential singular-element dis- 
tribution that is to be used in a certain case. 

As stated before, the purpose of this work was 
to identify (even if only qualitatively) a possible 
mechanism that could induce rolling moments on WAF 
configurations at zero angles of attack. It was 
therefore decided to use the simplest and least 
expensive available numerical approach. Conse-^ 
quently, the subsonic flow field over WAF config- 
urations was solved at zero angle of attack only. 

'I’li is restriction greatly simplxfied the modeling of 
the body and saved computation time. Xt did not 
severely limit the applicability of the calculated 
results since most WAF missiles fly at rather small 
angles of attack and the rolling moments at these 
angles are practically constant. 8 It was assumed, 
and later justified by the computed results, that 
at zero angle of attack the solution could be 
uncoupled. Since the fins were not generating any 
lift (except for the minute amount necessary for 
the induced rolling moment), their effect on the 
flow field, and consequently on the solution for 
the body, is believed to be negligible. It was 
therefore decided to simulate the body only first 
by a distribution of sources along its axis, and 
then to simulate the fins in the presence of the 
body by a modification of die nonlinear vortex- 
lattice method of Ref. 11. This uncoupled solution 
minimized the computation time, because the body- 
induced flow field had to be computed only once, 
and then it could be used for all subsequent fin 
arrangements . 

Body Modeling 

A good approximation of the body- induced flow 
field could have been obtained using one of the 
now-classic source-panel methods (such as discussed 
in Ref. VI). An even better approximation could 
have boon obtained by the simultaneous solution for 
the complete interaction flow field of the body and 
the wraparound surfaces with the method of Ref. 13* 
This method also would have accounted for the 
effects of the fins and their wakes on the body at 
finite angles of attack. However, such methods 
are time-consuming and unjustified for a qualita- 
tive investigation such ns this and, as later 
observed, not necessary at zero angle of attack. 

In this work the actual description of the 
ax isy mine trie body at zero incidence is by a dis- 
tribution of a finite number of discrete point 
sources along its axis* The intensities of those 
sources are determined by satisfying the flow- 
t a ago ney bounds ry condition [ E q . ( 2 ) ] on the same 
number of collocation points on the body. The 
details of this well-established method are given 
by von Kdriinln 1 ** and are only .summarized here. 
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Let (Q(I)} be the vector of the source inten- 
sities, and let fx c (I)} designate the locations 
along the body axis of these sources. The velocity 
at any given point (x,r) in the cylindrical system 
(x,r,G), induced by the sources system, is given by 


x - x c (I) 


u x (x,r) - E QW {[x _ Xc(I) p + 


u r (x,r) * 2 Q(I) 


I?L {[x - x c (X) 1 2 + r 2 


► 

77/2 

J J 


(5) 


where M Q is the total number of sources required 
to describe the body. The flow- tan gency condition 
[Eq, (2)] on the body surface is 


dr H B 11 r 

dx U OT + u x 


( 6 ) 


Equations (5) and (6) combine at any one of the 
collocation points (xp(J)} on the body surface 
into the following matrix equation: 

rr B (J) - [dr B (J)/d*ltx p (J) -* e (I) ll 

[ {[*p(J) - s c (I)) a + r B (J) J ) 3 /2 J Q '' 

** v'“ ir (J) } M s <?> 

Satisfaction of the boundary condition [Eq. (7)3 at 
one of the control points may be replaced by an 
imposed stagnation condition at the body apex to 
improve the accuracy of the solution in this 
region. 15 The linear system Eq. (7) is solved by 
the Gauss elimination method for the source inten- 
sities {Q(I)}. With the intensities of the sources 
known, the velocities induced by the body can bo 
calculated nt any point in the flow field using 
Eqs. (5). 

Obviously, the description of the body by a 
source distribution along Its axis is valid for an 
axisymmetric body at zero incidence only* For any 
other case, a more sophisticated scheme (such as 
that of Ref* 12) must be used; however, this 
description sufficed for the needs of this inves- 
tigation. Furthermore, the source distribution 
was computed only once in the absence of the fins 
(in order to save computation time) by assuming 
that the interaction of the fins with the body was 
negligible. This interaction was evaluated sep- 
arately from the results of the complete program, 
and the assumption was verified when it was found 
that the velocities induced by the fins on the 
body at zero incidence were negligible. This was 
expected because, at this condition, the fins do 
not generate lift and therefore do not disturb 
the flow. 

Wake Modeling 

As mentioned before, when simulation of wrap- 
around surfaces in uniform flow failed to produce 
rolling moments, the authors deduced from the roll 
reversal at sonic speeds and from the opposite 
sense of the rolling moment when the WAF-missile ’ s 
rocket motor was firing, that the wake structure 
must have nn upstream influence (in subsonic Clow) 
on the fins. A radial velocity component induced 
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on the antisymmetric fins could generate an asym- 
metric load on them and, therefore, would also gen- 
erate a rolling moment. 

With this in mind, it became necessary to simu- 
late not only the body, but also its wake. Since 
the potential flow model used could not account for 
a wake (neither the viscous wake of the body nor the 
plume of a rocket motor), the problem could be 
evaded by substituting a solid body for the inner 
wake enclosed by the dividing streamline. -For this 
to be done, the shape of the dividing streamline 
had to be known from external sources, and the 
resulting "solid body" was incorporated into the 
body simulation as if it were its inherent rear 
part . 


Three types of wake shapes are simulated in 
this work: 1) the free-undisturbed wake of a body 

in free flight; 2) the wake of a sting-mounted body 
in a wind tunnel; and 3) the plume of an under- 
expanded jet firing from a rocket-motor exit nozzle. 
Wake shapes are not readily available in any single 
reference, and a comprehensive literature review 
had to be undertaken in order to compile the neces- 
sary data. 2 * * * * * * * * * * * * * 16 In view of the length of the present 
paper, only the relevant conclusions of this review 
are given. 

1) The dividing streamline of a free wake was 
modeled by an ogival solid body. This shape was 
based on the experimental observations in' Refs. 17 
through 20, The flow visualization discussed in 
Refs. 17, 18, and 19 showed the separating stream- 
line to leave the body tangentially to the surface, 
and to enter the rear stagnation point at an angle 
between 30° and 40°. The tangent ogive used in the 
present model offered a good simulation of this 
pattern (Fig. 2). The most important characteris- 
tic of the wake — namely, its length to the rear 
stagnation point and its dependence on the Mach 
number (Fig. 2) — was based on Ref. 18. The 
Reynolds number did not affect the length and shape 
of the wake once the wake was turbulent . 18 1 2 0 The 
results of the experiments in Ref. 20 with differ- 
ent base shapes indicated that the presence of the 
fins should have no effect on the shape and total 
length of the wake. When the computer code was 
run, the influence of variations in the wake shape 
were investigated. Various wake shapes (e.g., 
parabolic or cubic polynomial) of equal length were 
tried, but the effects of these variations on the 
results were insignificant. 16 

2) Verification of the proposed mathematical 

model required a comparison with experimental data. 

However, such data were obtained from sting- 

mounted, wind-tunnel models. The sting changed 

the wake into that of a backward-facing step, vary- 

ing from an axisymmetric step 19 * (large step-height- 

to-body-diameter ratio) to an almost two-dimensional 

step 17 * (when the step is small). The location of 

the rear stagnation point (or reattachment point) 

on the sting has not been previously reported. The 

authors determined its approximate position from 

measured pressure distributions 21 and flow visuali- 

zations, 13 and defined an equivalent ogival after- 

body by an extrapolated rear stagnation point on 

the axis (Fig. 2) . The dependence of the length 

of this equivalent wake on the flow Mach number is 

shown in Fig. 2 for a body-to-sting-diameter ratio 

of 3:1. This ratio was chosen because it is widely 

used by wind-tunnel operators, and because it was 

used in the Arnold Engineering and Development 


Center (AEDC) 4-Ft Wind Tunnel subsonic and tran- 
sonic experiments reported by Dahlke. 9 Although 
the method of determining the position of the rear 
stagnation point is fairly crude, the influence of 
its inaccuracy is not large. The change in the 
induced rolling moment caused by a 20% variation 
in the length of the wake was in the range of 1-2% * 

3) The third wake simulated was a rocket- 
exhaust jet plume. WAF configurations are rocket- 
powered, and the induced rolling moments during 
the powered phase of the flight are also of inter- 
est. 6 The literature review of Ref. 16 defined a 
large number of parameters that determine the 
shape of the plume. Since the prime objective of 
this work was to obtain qualitative data on the 
induced rolling moments, and since experimental 
data on these moments during powered flight were 
not available for comparison, the plume was modeled 
for a limited number of parameters only. These 
parameters were for a sonic jet, blowing from a 
cylindrical nozzle exit into a subsonic outer flow. 
The jets were underexpanded, with jet-to-free-stream 
pressure ratios from 1 through 6. The typical geo- 
metrical features of the jet plume (the length of 
the first jet "diamond," its upstream vertex angle, 
and the diameter of its downstream end) under these 
conditions were taken from Refs. 22 and 23. With 
these three parameters specified, the shape of the 
solid body simulating the jet plume was modeled by 
a second-order polynomial. 16 


Fin Modeling 


With the body-induced flow field known, the 
pressure loading on the fins can be calculated and 
integrated to produce the induced rolling moment. 

In the present work, this is done by a nonlinear 
vortex-lattice method (VLM). 11 Each fin is divided 
into quadrilateral panels (Fig. 3), and bound 
horseshoe vortices of unknown strength are dis- 
tributed at the quarter chord positions of the 
panels (inset A in Fig, 3). The trailing vortices 
from the system are free to shed from the fins’ 
trailing and side edges, and after determining 
their equilibrium position in space, they form the 
vortex wake of the fins. The intensities of the 
vortices are again determined by the tangency- 
boundary condition [Eq. (2)]. This condition is 
imposed at a finite number of control points 
located at the midspan of each panel’s 3/4 chord 
lint. This tangency condition is satisfied by the 
mutual cancellation of the velocity components due 
to the free stream, body, and vortex flow fields 
that are normal to the fins’ surfaces at the con- 
trol points. 


4* -4 

The velocity dV, induced by a segment d£ 
of a line vortex of strength y at position r, 
is given by the Biot-Savart Law: 


^ _ Y dt x r 
4^ r 3 


( 8 ) 


When the contributions of all the vortices are 
integrated, the total induced velocity at any given 
point can be calculated. The component V n of the 
total induced velocity that is normal to the fin’s 
surface at a set of control points located at each 
panel’s 3/4 chord is given by 

Nv 

v n (J) = £ «(J.I)Y(X) j = 1 , 2 , . . .,N V ( 9 ) 
1=1 
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where N v is the total number of vortex panels 
with strengths y(I), and the matrix [H] is the 
influence coefficient matrix of all vortices (1) 
on points (J) . An influence coefficient is defined 
as the component, normal to panel (J) at its con- 
trol point, of the velocity induced there by a 
vortex of unit strength at (I) . The tangency 
boundary condition at the control points is satis- 
fied when the velocity component normal to the 
surface (that is, induced by the vortices) is can- 
celled by the normal component of the combined 
uniform and body-induced flow fields 

{U n } + [H] {y} = 0 (10) 

Satisfying this tangency boundary condition 
[Eq. (10)] at a given number of control points 
results in the same number of linear algebraic equa- 
tions for the vector of the unknown intensities of 
the vortices {y}. The solution of this system of 
equations requires, however, prior knowledge of the 
trajectories of the free vortices in order to 
account for the velocities induced by them on the 
control points. The trajectories, in themselves 
part of the solution, are either assumed to be known 
in the linear VLM codes or computed by an iterative 
process in the nonlinear VLM codes. 13 In the latter 
case, an initial guess for the trajectories is 
needed. It is usually assumed that the free vor- 
tices trail off to infinity as straight lines. 

With the intensities of the vortices known, the 
free vortices are allowed to follow local stream- 
lines. The vortex intensities are then recomputed 
and the process is reiterated until the vortex wake 
relaxes to its equilibrium rolled-up position. 

This nonlinear procedure was also followed in the 
present investigation. 16 However, it was found 
that since the fins were at zero incidence and the 
normal forces generated by them were of second 
order, their influence on the wake roll-up was 
almost insignificant. Consequently, the contribu- 
tion of the nonlinear process to the induced rolling 
moments was negligible, and they could be computed 
by the linear version of the code (its first itera- 
tion only) at a considerable reduction (75%) in 
computation time. 

With the strengths of the vortices known, the 
pressure distribution acting on the fins is com- 
puted by the Kutta-Joukowski theorem and is inte- 
grated relative to the body axis to give the rolling 
moment. 


Computational Results 

The main purpose of this work was to identify 
a possible mechanism inducing rolling moments on 
WAF configurations at zero angle of attack. The 
preliminary study indicated that, while no such 
moments could be produced in a uniform flow field, 
the radial velocity component at the base of the 
configuration (inward over a closed wake or outward 
over a jet plume) could induce rolling moments by 
interacting with the antisymmetric fins. 

To verify this hypothesis (even if only quali- 
tatively), the computation results had to be com- 
pared with experimental data. No such data were 
available for WAF configurations in free flight so 
wind-tunnel data generated by the TTCP (Ref. 7) had 
to be used. This was the reason for the simulation 
of the wake of a sting-mounted model. There were, 
however, drawbacks to the use Of these data. Not 


having been aware of the possibly decisive role of 
the wake flow in this case, the participants of 
TTCP did not specify in their reports (e.g.. 

Refs. 8-10) the model- to- sting-diameter ratios. 

One cannot be certain that the same ratio was used 
throughout the program. As shown in Fig. 2 this 
ratio has a first-order influence on the length of 
the wake and consequently also should have a non- 
negligible effect on the radial component of the 
velocity and on the rolling moments it induces. 
Since the experiments reported in Ref. 9 were done 
with a body-to-sting-diameter ratio of 3:1 (which 
is the common practice in wind-tunnel work) , and 
because of lack of information to the contrary, it 
was assumed here that all the data were obtained 
using this ratio. This assumption may have con- 
tributed to the scatter in the data, which already 
was considerable because of the difficulty in 
accurately measuring the relatively small rolling 
moment, especially at the lower subsonic Mach 
numbers. 9 The low accuracy of the experimental 
data limited the comparison with the computational 
results more to qualitative evaluation than to full 
verification of the mathematical model. 

The Standard Model 

The TTCP standard WAF configuration (Fig. 1) 
designated as B1F1 in the literature 9 was modeled 
for the comparison. The configuration was composed 
of an axisymmetric body (a 2-diam-long ogive nose 
and an eight-caliber cylindrical body) and a cru- 
ciform of fully deployed wraparound fins of rec- 
tangular planform that completely encircle the body 
when folded (Figs. 1 and 3). 

A distribution of a minimum of 50 sources was 
necessary for a reasonably smooth modeling of the 
body and its wake. The accuracy of the body model- 
ing was controlled by a comparison of the resulting 
body radius (of the zero streamline) between con- 
trol points with the actual radius, and by a com- 
parison of the local computed flow direction at 
these points with the body slope. The number of 
the sources and their spacing along the axis were 
automatically adjusted until prescribed accuracy 
criteria were satisfied. 

The modeling of the fins was done with 
100 uniform vortex panels (Fig. 3) per fin. This 
fine division guaranteed the numerical accuracy of 
the results. For instance, an increase of the 
number of panels by 30% changed the computed roll- 
ing moment by less than 0.5%. The computed pres- 
sure distribution over a fin that is shown in 
Fig. 4 indicates that the number of panels could 
be decreased and computation time could be saved 
by a nonuniform spacing of the panels with a far 
lower density on the upstream three-quarters of 
the fin and a denser division on its aft quarter 
where the pressure gradients are large. Another 
method to save computation time at a penalty of a 
2-3% error in the roiling moment was to model only 
one pair of fins on the configuration and multiply 
the computed moment by two. A linear solution for 
one pair of fins took about 125 sec on an 
IBM 370/165 computer; a cruciform fin configura- 
tion required about 500 sec. 

The linear solution needed a prescribed shape 
of the fins' vortex wake. Inset B in Fig. 3 
describes the options for this shape. The first 
assumption is that the free vortices lie on a 
truncated conical surface trailing the fin. The 
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user can specify the length and vertex angle of 
this cone and even a spiral twist of the vortices 
over the cone as a simulation of the wake roll-up. 

A nonlinear computation of the full-wake roll-up 
(found to be unnecessary for the present purpose) 
would require considerably more computer time. 

Comparison with Experimental Data 

In Fig. 5 the variation of the rolling-moment 
coefficient for the standard WAF configuration with 
the Mach number is compared with the experimental 
data of Ref. 10. These data were obtained with the 
fins in a cruciform (+) position and in an (x) posi- 
tion, and also with special miniature fin balances 
in addition to the main integral balance. The 
agreement between the computed results and the 
experimental data is good at Mach numbers 0.5 
and 0.8. The discrepancy at M = 0.3 cannot be 
explained by the present model. The general trend 
of the declining rolling moment with increasing 
Mach number is correctly obtained. The Goethert 
compressibility correction cannot be carried past 
M = 0.8> and the curve is extrapolated (dashed 
line) to the vanishing value measured at M = 1.0. 
Although not computable, this vanishing of the 
rolling moment at sonic speeds can be qualitatively 
explained by the present mathematical model. As 
the flow over the base of the body turns super- 
sonic, the turning of the flow (either inward or 
outward in the wake region) is centered on the base 
corner, and cannot exert any influence upstream on 
the fins. Thus, the fins experience only parallel 
flow and, consequently, no rolling moment is 
developed. Comparisons with wind-tunnel results 
of other WAF configurations are given in Ref. 16 
and show the same kind of agreement as is shown in 
Fig. 5. Also shown in Fig. 5 for comparison pur- 
poses is the rolling-moment coefficient for the 
same configuration in free flight. The shorter and 
more rapidly closing wake in free flight almost 
doubles the computed rolling moments. These 
results could not be verified experimentally for 
lack of such data. 

Other comparisons that were possible with 
experimental data are shown in Figs. 6 through 8. 
Figure 6 compares the rolling-moment coefficients 
at M a 0.8 for the standard configuration with 
the fins at different opening angles (Bq), where at 
0 D - 120° the fins are completely folded, and at 
0 D -0° they are fully deployed. Again, good 
agreement is observed, except for opening angles 
of 25° and 45°. No explanation can be offered for 
the discrepancy at these two angles. The same code 
was also used to compute the rolling moments gen- 
erated by a small antisymmetrical deflection of die 
fins at an angle (6) . The results are linear 
(Fig. 7). The computed slope of the curve, 

(C £ ) £ - 0.131/deg, falls close to the experimental 

value of (Cji)^ = 0.138/ deg. The final comparison 
was of the computed roll-damping coefficient (C$ ) 

with experimental data. To compute this parameter, 
the whole flow field was spun around the configura- 
tion with a velocity distribution adopted from 
two concentric cylinders with the outer cylinder 
at infinity. The results are compared in Fig. 8 
with data from the AEDC and Jet Propulsion Labora- 
tory (JPL) experiments reported in Ref. 9. The 
dashed line describes the results of the present 
mathematical model, and the solid line is obtained 
from these results by adding the roll damping due 


to the body boundary layer. The results are in 
fair agreement with the AEDC data (circles) . Com- 
parison with the JPL data is inconclusive because 
of the large uncertainty levels reported for these 
data (the AEDC data had no uncertainty levels 
shown) . 9 

Parametric Study 

Several other parameters were investigated, 
with details shown in Ref. 16. Some examples 
follow. 

Since the main assumption was that the base 
flow was responsible for the induced rolling 
moments, moving the fins upstream should weaken 
this moment as was demonstrated in Fig. 9. Roll- 
ing moments for fins flush with the base plane are 
compared with those for fins moved upstream half a 
body diameter and one body diameter, respectively. 
The rolling moment is clearly decreasing with 
increasing distance from the base, and extrapola- 
tion seems to indicate that it might vanish at 
some upstream distance of one and one-half to 
two body diameters from the base. It is interest- 
ing to note that the wraparound fins of Vought's 
MLRS missile are positioned upstream of the 
base 24 ; however, no explanation was given for this 
design. 

The strong effect of the proximity to the 
base is demonstrated also in Fig. 10, where the 
influence of fin leading-edge and trailing-edge 
sweep angle was studied. While leading edge 
sweepback reduces the rolling moment only slightly, 
an identical forward sweep of the trailing edge 
sharply reduces the rolling moment. This result 
could be anticipated because most of the flow 
turning induced by the base flow would be felt by 
the aft portion of the fin as was indeed demon- 
strated by the pressure distribution in Fig. 4. 
Moreover, most of the load there is carried by the 
outer aft corner of the fin that is sliced off 
when the trailing edge is swept forward. While 
there are no experimental data available for 
swept-forward trailing edges, Ref. 9 reports 
results for swept-back leading edges. However, 
the extremely large scatter in these data prevents 
any significant comparison. A possible experimental 
indication of the influence of the wake can be 
found in the strong effects that changes in the 
geometry of the fin trailing edge had on the roll- 
ing moment , 2 b 

The effect of the span of the fins is shown 
in Fig. 11. In certain missile designs the span 
of the fins is increased by increasing the sector 
angles (HO covered by the fins and by making them 
overlap (an arc angle of 90° is the common fin 
full span) . As expected, the larger area of the 
fin generates a higher rolling moment, but the 
increase is not linear. 

A further demonstration of the strong effect 
of the wake flow on the fins is presented in the 
rolling moments due to a jet plume of pressure 
ratio Pj/Pco “ 6 (Fig. 12). The results are shown 
for the only case of jet plume that was simulated 
(see the section on wake simulation) at the highest 
pressure ratio, because the rolling moment is then 
at its peak. For lower pressure ratios, the 
moment decreases in size and vanishes for an 
ideally expanded jet. As shown in Fig. 12, the 
outward radial component of the wake flow changes 
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the direction of the rolling moment, which is now 4 Sawyer, W., Monta, W., Carter, W., and 


about four times as great as the moment induced by 
the inward flow of the closed wake. 

For additional examples (e.g., canard WAF, 
other asymmetric fin configurations, swept fins, 
etc.) see Ref. 16. 


Conclusions 

By using a source-distribution simulation for 
the body and a vortex-lattice simulation for the 
wraparound fins, it is shown that the radial veloc- 
ity component induced on the fins by the wake flow 
field could be responsible for the rolling moments 
encountered by WAF configurations in subsonic flow, 
at zero angle of attack. No other mechanisms were 
found; however, their existence cannot be ruled out. 
The present simulation estimates rolling moments of 
the right magnitude and correct trends. 

Assuming that the mechanism described by the 
present mathematical model is valid, then the 
induced rolling moments could be eliminated by 
mounting the fins about two body diameters upstream 
of the base of the body. Another method to prevent 
roll reversal in the accelerating phase is to 
accelerate into supersonic speed before the rocket 
motor burns out, because the roll direction is the 
same in powered subsonic flight and in supersonic 
flight . 3 This, however, will not solve the problem 
when the WAF configuration slows down to subsonic 
speeds at the end of its flight. The negative roll- 
ing moment experienced by WAF configurations in 
supersonic flow cannot be explained by the present 
purely subsonic mathematical model. It may be due 
to shock wave focusing on the concave side of the 
curved fin with the resultant pressure difference 
acting in the direction of the convex surface. 

This hypothesis requires further theoretical work 
as does a parametric study of various effects on 
the rolling moments. The authors also intend to 
apply the source-panel method 13 to describe the 
body and the nonlinear vortex-lattice method for 
the fins in order to compute WAF aerodynamics at 
higher angles of attack. 
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Fig. 3 Schematic modeling of WAF configurations. 16 




Fig. 2 Modeling of the closed wake in free flight Fig. 4 Load distribution over fin on standard WAF 
and on a sting mount. configuration at M = 0 and zero angle of attack. 
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Fig. 5 Rolling- moment coefficient, compared with Fig. 8 Roll-damping-moment coefficient, including 

experimental data (standard WAF configuration). body contribution, compared with the experimental 

data. 



Fig. 6 Rolling-moment coefficient vs. fin-opening Fig. 9 Free-flight rolling-moment coefficient 
angle compared with experimental data. variation with fin axial position. 



Fig. 7 Rolling-moment coefficients due to fin- 
differential deflection. 
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Fig. 10 Free-flight rolling-moment coefficient 
variation with leading-edge and trailing-edge sweep 
angle. 




Fig, 11 Free- flight rolling-moment coefficient vs. 
fin span angle. 


Fig. 12 Rolling-moment coefficient with a 6:1 
pressure ratio jet plume. 
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